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Abstract 

We propose a numerical criterion which can be used to obtain accurate and reliable 
values of the ordering temperatures and critical exponents of spin glasses. Using this 
method we find a value of the ordering temperature for the ± J Ising spin glass in three 
dimensions which is definitely non-zero and in good agreement with previous estimates. 
We show that the critical exponents of three dimensional Ising spin glasses do not appear 
to obey the usual universality rules. 

The full explanation of the universality rules for critical exponents in second order transi- 
tions through the renormalization group theory is one of the most impressive achievements 
of statistical physics. The universality rules for such transitions state that the critical ex- 
ponents depend only on the space dimension d and a few basic parameters : the number 
of order parameter components n, the symmetry and the range of the Hamiltonian |]]. No 
other parameters are pertinent. In fact it is known that there are exceptions to universality - 
in certain two dimensional (2d) Ising systems with regular frustration the critical exponents 
vary continuously with the value of a control parameter 0]. As far as we are aware, no 
results of this type have been reported in any three dimensional (3d) family of Ising systems; 
it has been tacitly assumed that non-universality is very exceptional. 

As compared to standard second order transitions, the situation concerning Ising Spin 
Glasses (ISGs) is much less clear; indeed the history of ISG simulations has been plagued 
by technical difficulties associated with long relaxation times. For two decades the very 
existence of a finite temperature transition in the 3d ISG has been hotly contested; as it 
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is obviously essential to have a reliable value of the ordering temperature before obtaining 
accurate critical exponent estimates, it has been difficult to make stringent numerical tests 
of universality in 3d ISGs. 

We will present a numerical criterion which can in favourable cases provide precise and 
reliable values for the ordering temperature T g and for the critical exponents of a spin 
glass, with a moderate level of computational effort. If an independent estimate of the 
ordering temperature is available the criterion leads to a convenient method for estimating 
the exponents. We study 3d ISGs with various sets of interactions and we conclude from the 
data that the 3d ± J interaction ISG has a well defined non-zero Tg which can be estimated 
accurately, and that universality in the usual sense does not hold in 3d ISGs. 

It would appear probable that glassy transitions in general have a much richer critical 
behaviour than have conventional second order transitions. 

Thus, technically the most difficult problem in numerical ISG studies is the correct iden- 
tification of the transition temperature T g . For the 3d ISG with random ±J near neighbour 
interactions on a simple cubic lattice, which has been the subject of a considerable ammount 
of work, T g has been estimated in two ways. Ogielski || studied in massive simulations the 
divergence of the spin glass susceptibility, of the correlation length, and of the relaxation 
time of the autocorrelation function 



q(t) = [< Si{t)Si(0) >] (1) 

in order to estimate T g and the critical exponents. However his analysis has been questioned 
because of the possibility of ambiguities in the manner of identifying a divergence, if non- 
conventional temperature dependencies are invoked |§] . Bhatt and Young f5j used a finite size 
scaling technique; they measured the Binder cumulant for the fluctuations of the equilibrium 
autocorrelation function 
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as a function of sample size L. The curves gi{T) for different L should all intersect at T g \ in 
the 3d ±J ISG case the curves indeed intersected but did not appear to fan out below the 
apparent T g . Only recently have intensive numerical studies shown that a weak fanning out 
at low temperatures really does occur || [7|]. Even with results of high statistical accuracy 
to hand, Kawashima and Young || give a number of caveats concerning the interpretation 
of their own data. 

We will describe an alternative criterion for determining T g . First, scaling rules tell us || 
that for a large sample in thermal equilibrium at T g the relaxation of the autocorrelation 
function takes the form 

q(t) = \t~* (3) 
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with the exponent x related to the standard static and dynamic exponents 77 and z through 



« = (4) 

Secondly, the out of equilibrium relaxation of two randomly chosen replicas A and B of 
the same sample towards equilibrium at T g depends on another combination of the same 
exponents ||. The out of equilibrium spin glass susceptibility is defined as 

X'soit) = [< S*(t)S*(t) > 2 ] (5) 

and it increases with time as 

t h with h = — —. (6) 

z 

Suppose we take {T}, a series of trial values for T g ; from measurements of x and h on large 
samples at each T we can deduce from equations f| and |^ a set of apparent or effective 
exponents 

z(T) = (8) 
K J 2x + h K ' 

Finally in another set of simulations on the same system at different [small] sample sizes L, 
from standard finite size scaling rules || for the fluctuations in the autocorrelation function 
in equilibrium at T g we have 

L d - 2 < q 2 > oc L~ v (9) 

If we again take a series of trial values of T g and fit the results using this form at each Tj 
we will obtain a second series of apparent exponent values ^(T). (This type of fit will only 
be appropriate close to and below T g ; at higher T another factor appears on the right hand 
side §). 

We now plot r/i(T) and 772 CO against T; for consistency the curves must intersect at the 
point [77, Tg] which represents the true critical exponent r\ and ordering temperature T g of 
the system. At this temperature and this temperature only the functional forms of equations 
[3], ^| and ^ should be exact; at neighbouring temperatures these forms are only approximate 
but close to T g they will be adequate to parametrise the numerical data. Once T g is fixed by 
the intersection we can obtain z using the z(T) curve given above, and with known r\ and T g 
we can go on to fit < q 2 > data for temperatures above T g to obtain the exponent v. From 
scaling relations, once we dispose of 77 and v all other static exponents are determined. 

We show in figure [l] estimates for r\\ (T) and 772 (T) for the 3d ±J ISG calculated using data 
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Figure 1: % (o) and ry 2 (•) vs T for various distributions, a) ±J, b) Uniform c) Gaussian 
and d) decreasing exponential. Note that the scale on the x axis is different for each plot. 
Error bars on individual r) points are about ± 0.02. 



taken from the literature : x(T) from ||, h{T) from || |§, and the spin glass susceptibilities 
for different assumed values of T g ; (T g = 1.0 from the data given in ||, T g = 1.11 from ||, 
and T g = 1.175 from ||). There is a well defined crossing point with T g = 1.165 ± 0.01 and 
7] = —0.245 ± 0.02. Using the curve for z(T) from equation [5] we estimate z = 6.0 ± 0.2. 

The values obtained in this way are at least as precise as previous estimates and are very- 
close to the central values given by Ogielski @ (T g = 1.175 ± 0.025, r] = -0.22 ± 0.05, 
z = 6.0 ± 0.8), corroborating his analysis. On the other hand the T g is marginally outside 
the error bars quoted by Kawashima and Young (T g = 1.11 ±0.04) who use extensive Binder 
cumulant data ||. The difficulty in applying this latter method to the 3d ± J ISG case is that 
the gi{T) curves lie very close together below T g so the intersection point is sensitive to small 
changes in individual Ql curves. Even with extreme statistical accuracy, small corrections 
to finite size scaling (invoked as a possibility in 0) can change the apparent position of the 
intersection point significantly. The results of ref || could be rendered consistent with the 
present analysis if the gi values for the smallest samples studied were affected by corrections 
to finite size scaling at the 1% level. 

The present method is much less sensitive to problems of systematics than are either of the 
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other techniques outlined above. First, both x and h are determined using "large" samples so 
finite size corrections should be unimportant 0, |9| . Secondly h is measured out of equilibrium 
and so is not subject to the problems of long equilibration times. The fact that no preparatory 
anneal is required also means that the measurements are economical in computer time. The 
measurements of x need careful equilibration but systematic tests using successsively longer 
preliminary anneals allow one to obtain reliable values. Numerical data []3], [M]] show that in 
ISGs q(t) already takes on the asymptotic form, equation [| from quite early times t ~ 2 
MCS (Monte Carlo Steps), and that sample to sample variations in the values of x are small 
so extensive averaging over very large numbers of samples (an essential condition for good 
gi data) is unnecessary. Thus the curve r]i(T) can be established accurately with moderate 
numerical effort and minimal systematic error. For the finite size scaling data from which 
772 (T) is deduced, thorough equilibration is necessary but by studying pairs of replicas |5j and 
again testing with increasing anneal times it is easier to obtain accurate values of < q(t) 2 > 
than the combination of moments which constitute the Binder cumulant. Again, the sample 
to sample variability is much less for < q(t) 2 > than for the Binder cumulant. In the 3d 
±J ISG the two curves r\\ (T) and 772(2^) intersect cleanly, figure [l], so the determination 
of the crossing point should not be very sensitive to minor deviations from scaling or small 
statistical uncertainties. Finally, no hypothesis is made concerning the way divergences occur 
except the essential assumption that standard scaling rules (as opposed to universality rules) 
hold. The excellent overall agreement between Ogielski's estimates H and the present ones 
gives considerable confidence in the general coherence of the standard scaling approach and 
appear to make any exotic scaling assumption unnecessary. 

We therefore consider that both f]i(T) curves can be calculated with little in the way of 
disguised systematic errors; as they stand the T g and exponent values that we quote should 
not only be precise but reliable. 

We have made further simulations on another 3d ISG with ± J interactions; this is the 



the fully frustrated system with 20% random bond disorder that we studied in | 10fl . We 
already established an accurate value of T g (T g = 0.96) for this spin glass from Binder 
cumulant measurements, and we now have measured the exponents x and h at T g together 
with an estimate of i] from the spin glass susceptibility (see Table [I]). The data are very 
consistent with each other and lead to an 77 value which is less negative and a z value which 
is smaller as compared with those of the standard ± J ISG. This difference already indicates 
the non-universality of these two exponents in 3d ISGs. 

We have also carried out extensive simulations on 3d ISG systems with different sets of 
near neighbour interactions. For the 3d ISGs with near neighbour Uniform, Gaussian and 



decreasing Exponential interactions (see |TlJ for the definitions of the distributions with the 
correct normalizations), the data are shown in figure [I]. Simulations were done on samples 
with L — 16 for x, L — 10 for h, and samples from L = 2 to L = 6 for < q(t) 2 >. 
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Careful anneals were carried out where appropriate, checked by the prescription given in 
0. At each temperature, 10 samples were used for x, 500 for h and 2000 to 200 depending 
on L for < q(t) 2 >. We estimate that the r)i(T) curves are on large enough samples for 
there to be virtually no finite size correction, so the values can be taken as definitive (apart 
from statistical errors), but measurements on larger samples could modify the i]2(T) curves 
marginally. It can be seen that the r)(T) curves again cross cleanly for the Uniform case with 
a more negative r\ than for the ± J case. However for the Gaussian and Exponential cases 
it turns out that the two curves are much more similar to each other making it difficult to 
identify T g precisely; for these distributions we have to fall back on an alternative method 
to estimate T g . 

The Migdal-Kadanoff (MK) scaling approach is known to give reasonable values of the 
ordering temperature for Ising spin glasses |T2|, [T3|, [TJj]. We have followed the particular 
method used by Curado and Meunier MM but with improved statistical accuracy. It turns 



out that with a scale factor 6 = 2 the MK estimate for the 3d ±J ISG T g is 1.16 ± 0.01, 
precisely the same as the value we have obtained above from the simulations. This perfect 
agreement is certainly fortuitous (though in 4d where the MK method should be much 
poorer, the disagreement in T g between the 6 = 2 MK estimate and an accurate simulation 
value is only 15% Ul5|), but we argue that as agreement happens to be excellent for the 
± J case, if we apply the same method with the same scale factor 6 to other 3d ISGs with 
different sets of interactions, we should obtain T g estimates which should again be very close 
to the real values. We obtain MK T g values which are 1.00, 0.88, and 0.72 for the Uniform, 



Gaussian and Exponential distributions respectively [To]. The Uniform distribution value 
is in good agreement with the simulation value and the other two T g values are within the 
range of T where the simulation curves for r/i(T) and ^{T) overlap. The Gaussian T g and rj 
are in good agreement with earlier estimates ||. Putting uncertainties at ±0.05 for possible 
systematic errors in the Gaussian and Exponential MK T g estimates, we obtain the set of 
exponent estimates shown in Table [3]. 

According to the usual universality rules, the form of the interaction distribution should 
not be a pertinent parameter as concerns the critical exponents. Here we find that 3d ISG 
systems which differ only by this distribution function show quite different r] and z values, 
Table [IJ. The results indicates a breakdown of conventional universality in 3d ISGs. 

In order to show that the apparent non-universality is not an artefact, we will turn back 
to the raw x and h data for the ±J and Uniform cases. In figure ^ we have plotted the 
values of these parameters as a function of T; the error bars are about ±0.005 for h and 
±0.002 for x. If universality holds 

h(TJU)) = h(T g (J)) (10) 
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Figure 2: h(T) and x(T) for ±J (A) and Uniform (o) distributions. The temperature scale 
is common. The dashed line corresponds to the example given in the text 
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System 


T 


x{T g ) 


KT g ) 


V 


z 


±J 


1.165 ±0.01 


0.064 


0.38 


-0.245 ±0.02 


6.0 ±0.2 


FFd0.2 


0.96 ±0.02 


0.091 


0.437 


-0.12 ±0.02 


4.85 ±0.3 


U 


1.05 ±0.03 


0.054 


0.41 


-0.375 ±0.03 


5.8 ±0.5 


G 


0.88 ±0.05 


0.035 


0.355 


-0.50 ±0.04 


7.1 ±0.6 


Exp 


0.72 ±0.05 


0.02 


0.275 


-0.62 ±0.12 


9.5 ±0.7 



Table 1: Temperature of transition and critical exponents for several distributions. The 
distributions are in order (i) random ±J interactions, (ii) Fully Frustrated lattice with 
20% disorder (Hi) random Uniformly distributed interactions, (iv) random Gaussian 
interactions, (v) random Decreasing Exponential interactions 



and 

x{T g {U))=x{T g {J)). (11) 

By inspection, whatever trial value T* we choose for T g (J) within the generous limits T* = 
1.0 to 1.3 provided by the figure, the relation [TO] leads to us to a T*{U) such that x(T*(U)) 



is considerably smaller than x(T*(J)). For instance with T*(J) = 1.16, T*{U) = 0.88, 
x{T*(J)) = 0.064, x(T*(U)) = 0.036. The data cannot satisfy and [II] simultaneously, 
demonstrating non-universality. 

For the 2d regularly frustrated systems which show continuous variation of critical expo- 
nents, the breakdown of universality is necessarily associated with the existence of a marginal 



operator [16] and it has been pointed out that when breakdown occurs, it does so in Ising sys- 



tems having more than two ground states |17] and hence with n, the number of components 



of the order parameter, greater than 1 [IE . On the Parisi image of finite dimension ISGs 
19|1, n is essentially infinite; it would be of interest to identify possible marginal operators. 
We can note that in the regularly frustrated 2d systems quoted above, v varies continuously 
but 7] is constant so "weak universality" p0| still holds. This is not the case for the randomly 
frustrated systems we have studied. 

It would appear that universality breakdown could be much more prevalent than was 
suspected, and it may well be the rule rather than the exception at spin glass or glass 
transitions. 
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